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ABSTRACT 

The BO. 2 V magnetic star rSco stands out from the larger population of massive 
OB stars due to its high X-ray activity, peculiar wind diagnostics and highly complex 
magnetic field. This paper presents the discovery of the first two r Sco analogues - 
HD 66665 and HD 63425, identified by the striking similarity of their UV spectra to 
that of T Sco. ESPaDOnS spectropolarimetric observations were secured by the Mag- 
netism in Massive Stars CFHT Large Program, in order to characterize the stellar and 
magnetic properties of these stars. CMFGEN modelling of optical ESPaDOnS spectra 
and archived lUE UV spectra showed that these stars have stellar parameters similar 
to those of T Sco. A magnetic field of similar surface strength is found on both stars, 
reinforcing the connection between the presence of a magnetic field and wind pecu- 
liarities. However, additional phase-resolved observations will be required in order to 
assess the potential complexity of the magnetic fields, and verify if the wind anomalies 
are linked to this property. 
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1 THE MAGNETIC B-TYPE STAR r SCO 

Very little is known about the magnetic fields of hot, massive 
OB stars. Due at least in part to the challenges of measure- 
ment, direct evidence of magnetic fields in massive stars is 
remarkably rare - there are perhaps 30 known examples - 
and few have been studied in detail. 

Even as a member of the elusive class of magnetic mas- 
sive stars, the BO. 2 V star r Sco is recognized to be a pe- 
culiar and outstanding object. First, r Sco is distinguished 
by its magnetic field, which is unique because it is struc- 
turally far more complex than the mostly-dipolar (multipole 
mode I = 1) fields usually observed in magnetic OB stars. 
According to the Z eeman Doppler Imaging performed by 
iDonati at all l|2006t ). the magnetic topology of r Sco presents 
significant power in multipole modes up to ? = 5 with a mean 
surface field strength of ~ 300 G, and the extrapolated 3D 



* Based on observations obtained at the Canada-France-Hawaii 
Telescope (CFHT) which is operated by the National Research 
Concil of Canada, the Institut National des Sciences de I'Univers 
of the Centre National de la Recherche Scientifique of France, and 
the University of Hawaii. 



field structure shows an intricate assortment of open field 
lines and closed magnet ic loops (illustrated in Fig. 11 and 
12 of lDonati et~alll2006l '). 

Secondly, r Sco also stands out from the crowd of early- 
B stars because of its stellar wind anomalies, as diagnosed 
through its odd UV spectrum and X-ray emission. The opti- 
cal and non-resonance UV lines of r Sco, that probe the con- 
ditions of the photosphere, are typical of an early-B dwarf 
star. However, the UV resonance lines arising from the stel- 
lar wind are not those typical of such a star. Figure [1] shows 
a comparison of the UV spectrum of r Sco (second spec- 
trum) with the spectrum of a normal BO main sequence 
star (bottom spectrum) and a BO giant (top spectrum). 
The resonance wind lines of Civ (right panel). Si IV and 
Nv (left panel) are indicated with dashed lines. The UV 
wind line morpholog y of normal, early B stars conforms to 
a 2-D spectral grid jWalborn et all 119951 ). Typically, Civ 
strengthens with increasing temperature and luminosity. N v 
is very weak on the main sequence at spectral type BOY 
but strengthens with temperature and luminosity. For stars 
with fixed Civ and Nv profiles, Siiv is strictly luminosity 
dependent and breaks the degeneracy, r Sco does not fit into 
this grid. It has strong Nv indicating that it should be well 
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Figure 2. Fits to tlie optical Hel-Il and H/3 lines for HD 66665 
(loft) and HD 63425 (right) 

above the main sequence. However, its C IV lines are only 
slightly stronger than typical and not distinctly wind-like, 
suggesting a near main sequence luminosity. Finally, its Si iv 
profiles are unique, a bit stronger than typical class V stars, 
but unlike normal, early giants. As a result, this stars lie 
outside of the normal classification grid, suggesting a more 
highly ionized outflow than typical. The hard X-ray emis- 
sion of r Sco also sugg ests hot plasma, in excess of 10 MK 
l|Cassinelli et al.lll994i ). We suspect that this is connected 
to the strange UV properties. Like its magnetic field, the 
stellar wind of r Sco possesses observational characteristics 
that make it unique among OB stars. 

Interestingly, the wind lines of r Sco have been shown 
to vary periodicall y with the star's 41 d rotation period 
l|Donati et al.ll2006l ). which is interpreted to mean that the 
wind is structured and modulated by the magnetic field. 
Clearly the magnetic field exerts an important infiuence on 
the wind dynamics. What is not clear is whether the wind- 
line anomalies described above are a consequence of the un- 
usual complexity of r Sco's magnetic field, a general conse- 
quence of wind confinement in this class of star, or perhaps 
even unrelated to the presence of a magnetic field. Because 
such wind anomalies have never been observed in any other 
star, magnetic or not, this issue has remained unresolved. 

The identification and analysis of additional stars with 
wind properties similar to r Sco would therefore represent 
an important step toward understanding the origin of these 
peculiarities. In this paper, we present two early B-type stars 
- HD 66665 and HD 63425 - that we identified to be the first 
r Sco analogues based on their UV spectra, which are strik- 
ingly similar to the UV spectrum of r Sco. Figure [T] shows 
lUE archives spectra of these two stars (third to fifth spectra 
from the top) . The discovery of wind anomalies naturally led 
to an investigation of the magnetic properties of these stars 
by the Magnetism in Massive Stars (MiMeS) collaboration 
l|Wade et al.ll201ll ). 



2 OBSERVATIONS 

Spectropolarimetric observations of HD 66665 and HD 63425 
were obtained with ESPaDOnS at the Canada-France- 
Hawaii Telescope. We acquired 4 high-resolution {R ^ 
65 000) broad-band (370-1050 nm) intensity (Stokes I) and 



Table 1. Journal of ESPaDOnS observation listing the UT date, 
the heliocentric Julian date (2 400 000-I-), the total exposure time 
and the s/n per 1.8kms~^ velocity bin at 540nm. Column 5 
gives the probability that the observed signal in Stokes V LSD 
profiles is due only to noise, and column 6 the global longitudinal 
field component obtained from the first moment of the Stokes V 
profiles. The available archival lUE spectra are also listed. 



Date 


HJD 


^cxp 


s/n 


FAP 




(UT) 




(s) 






(G) 






HD 66665 






2010/02/26 


55 253.97 


4 400 


728 


< 10-8 


-93 ±4 


2010/03/02 


55 257.98 


1200 


236 


< 10-8 


-3 ±15 


2010/03/04 


55 259.98 


1200 


330 


< 10-8 


9 ±10 


2010/03/05 


55 260.98 


1200 


276 


2 X 10-** 


41 ±13 


1983/12/04 


45 436.77 


900 




lUE swp 21680 


1987/03/11 


46 866.49 


420 




lUE swp 30484 


1987/03/12 


46 866.58 


1080 




lUE swp 30486 






HD 63425 






2010/02/28 


55 255.84 


1200 


520 


< 10-8 


117±8 


2010/03/03 


55 258.88 


1200 


289 


< 10-8 


113 ±15 


2010/03/07 


55 262.76 


1200 


435 


< 10-8 


132 ±10 


2010/03/08 


55 263.83 


1200 


361 


< 10-8 


125 ±12 


1987/03/12 


46 866.53 


600 




lUE swp 21679 


1983/12/04 


45 672.74 


720 




lUE swp 30485 


^ A detectable 


Stokes V 


signal 


can be 


observed 


with high- 



reso lution instr uments even if the global longitudinal field is null 
fsee lPetitllioTlh . 



Table 2. Summary of stellar properties of tSco, HD 66665 and 
HD 63425. 





tSco^ 


HD 66665 


HD 63425 


Spec, type 


B0.2V 


B0.5V 


B0.5V 


Teff (kK) 


31 ± 1 


28.5 ±1.0 


29.5 ± 1.0 


logff (cgs) 


4.0 ±0.1 


3.9 ±0.1 


4.0 ±0.1 


R* (Rq) 


5.6 ±0.8 


K c +3-3 
5.5 -2.7 


6.8 -2.0 


log(L./LQ) 


4.47 ±0.13 


+0.40 
4.25 -0.60 


. ^„ +0.46 
4.50 -0.30 


(Mq) 


11 ±4 


9 -7 


17 -9 


V sin i (km s— '^) 


< 13 


< 10 


^ 15 


M (lO-^Mgyr-l) 


61 _2 


< 0.45 


< 0.75 


Vao (kms-^) 


~ 2000 


~ 1400 


~ 1700 



Paramters from ISimon-Diaz et al ] l|2006h . and iMokiem et al.l 
1I2OO5I) for M and v^. 



circular polarisation (Stokes V) spectra for each star, dur- 
ing about 10 consecutive nights in Feb/Mar 2010. A sum- 
mary of the ESPaDOnS observations is given in Table [1] 
For a complete description of observation procedure and re- 
duction proce dure with the l ibre- esprit package provided 
by CFHT, see iDonati et al.l (|l997l ). In addition to the ES- 
PaDOnS observations, the available archival lUE spectra 
are also listed in Table [T] 



3 STELLAR PARAMETERS 

In order to determine the stellar and wind parameters of 
HD 66665 and HD 63425 we used non-LTE model atmo- 
spheres and spectrum synthesis from the CMFGEN code. CMF- 
GEN solves for the radiative transfer equation in a spherically 
symmetric outflow together with the radiative and statisti- 
cal equilibrium equations, taking into account effects such as 
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Figure 1. lUE spectra of rSco and its analogues HD 66665 and HD 63425 (second to fifth spectra). For comparison, typical spectra 
of a BO dwarf and a BO giant are shown (bottom and top respectively). The dashed lines indicate the wind lines NvAA1239, 1243, 
SilvAA1393, 1403 and ClvAA1548, 1550. The best-fit synthetic spectra, from which stellar parameters were derived from for HD 66665 
and HD 63425, are shown in red. 



line-b lanketing and X-rays (for details, see lHillier fc Milleij 
1 19981 ). 

We used optical spectra to derive the surface gravity 
(log(?) and effective temperature (Tcff), using classic line 
diagnostics (i.e., Balmer line wings and Hel-ll lines). UV 
spectra from lUE were used to infer stellar wind proper- 
ties. All results are summarized in Table [2l along with those 
of rSco for comparison. A solar chem ical composition was 
adopted, following Asplund et al.l (|2009iy We postpone a de- 
tailed analysis regarding chemical abundances for a future 
study. 

3.1 HD 66665 

The various estimates for the distance of HD 66665 range 
from about 1 to 2 kpc (Savcigg c t al. 1985; Savage & Massa 
Il987l : iDiplas fc Sawel 1 1993 : Ivan LeeuwenI feoOTi ). We 
adopted an intermediate value of 1.5 kpc and assumed an 
uncertainty of a factor of twcQ. The stellar luminosity 
was then derived by fitting the spectral energy distribution 
using the lUE spectrum and UBVJHK photometry, con- 
verted to flux points. An E{B - V) = 0.02 and R = 3. 1 
were used to fit the observed data jPiplas fc Savagelll994 ). 
The stellar radius (R*) and mass (M*) follow directly from 
the chosen L*, Tofi, and log p. 

A very good agreement is obtained for the optical di- 
agnostic lines in HD 66665. In Fig. [2] (left) we present fits 
to different Hei-ii lines and H/3. Our model was convolved 
with a V sin i of 8 km s~^ . We note however that this value is 
an estimate. Values in excess of lOkms"^ tend to produce 
profiles which are generally broader than the ones observed 
but accurate lower values are difficult to determine. No sig- 



^ Note that the recently revised paralla x for HD 66665 is v ery 
uncertain, namely, n = —0.94 it 1.05 mas jvan Leeuwen|[2007^ . 



nificant variability was observed in the optical spectra of 
HD 66665. 

The UV spectru m of H D 66665 was first pointed out as 
unusual by .Grigsbvl (|l992l ) . The absorption shown by C iv 
cannot be reproduced in detail by our models (see Figure 
[1]). Several tests were performed, but the observed profile 
is always too deep compared to the synthetic ones. On the 
other hand, we could achieve an acceptable fit for the Siiv, 
and N v lines. The inclusion of X-rays in our CMFGEN models 
was essential in producing the Nv feature, since the wind 
structure is not sufficiently ionized without this component. 

Despite the problem found for C IV, we could derive an 
upper limit for the mass-loss rate (M) from the emission part 
of the synthetic P-Cygni profile of this same line. When M is 
greater than about 4.5 x 10~^° M0 yr~^, the model predicts 
a relatively strong P-Cygni emission which is clearly absent 
in the observed spectrum, as it is in the spectrum of r Sco. 
Furthermore, only a rough estimate for the wind terminal 
velocity («txj ~ 1400kms~^) could be made. 

3.2 HD 63425 

We adopted a d istance to HD 6342 5 of 1136 pc (from n = 
0.88 ± 0.36 mas: Ivan LeeuwenI [20071 ). The stellar luminosity 
was then derived in the same way as in HD 66665. For the 
reddening parameters, w e use i? = 3.1 and E{B — V)= 0.10 
(|Savage fc Massal Il987l ). Our fits to the Hel-ll diagnostic 
lines and H/3 are presented in Figure [2] (right). The optical 
spectra of HD 63425 and HD 66665 are similar. However, our 
analysis indicates a slightly higher effective temperature and 
surface gravity for HD 63425. Moreover, Balmer line cores in 
HD 63425 are shallower, especially Ha and H/3. Line filling 
by a much higher mass-loss rate than the one in our final 
model is not a reasonable explanation here, since it would 
produce strong UV wind profiles which are not observed. 
A possible explanation for this line filling could be contam- 
ination (emission) from material trapped by the magnetic 
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field, distinct from the outflowing wind. This problem does 
not occur in HD 66665, where we have good fits to entire 
Balmer line profiles. 

The apparent rotational velocity of HD 63425 is low, 
but an accurate value is difficult to infer from the observed 
spectrum. Values above 15kms~^ tend to produce profiles 
that are broader than the ones observed. After various tests, 
we adopted a wsinz of llkms"^ for the convolution of our 
final model. 

There are only two lUE spectra of HD 63425. Inter- 
estingly, although the optical spectra of HD 63425 do not 
show any variation, these UV spectra (obtained over 3 years 
apart) reveal large variations in the wind lines (see the two 
spectra shown in Fig.[T]). The SWP 21679 spectrum presents 
these features larger and deeper than the ones found in 
SWP 30485. An interesting possibility was proposed to ex- 
plain a similar issue in the star r Sco. The passage of the 
magnetic equator across the line-of-sight would reveal ma- 
terial trapped by t he magnetic fi eld, providing the extra ab- 
sorption observed jPonati et al ][2006). Unfortunately, given 
the scarcity of UV data for HD 63425 we have no means to 
test this idea at present (for example, by verifying the exis- 
tence of absorption episodes) . 

We chose to focus on the SWP 30485 spectrum, where 
we have shallower absorption profiles. Nevertheless, the 
mass-loss rate obtained can be safely used as an upper limit. 
Our fit is illustrated in Figure [T] Again, the use of X-rays 
was mandatory in producing a Nv feature comparable to 
the observed one. We note that the C iv AA1548, 1550 ab- 
sorptions could be partially or entirely due to the interstel- 
lar medium. We summarize the wind parameters obtained 
in Table |2l which are similar to those of HD 66665. 



4 MAGNETIC FIELD 

In order to increase the magnetic sensitivity of our data, we 
applied the Le ast-Squares D e convo lution (LSD) procedure, 
as described bv lDonati et al.l l|l997h . This procedure enables 
the simultaneous use of many lines present in a spectrum 
to detect a magnetic field Stokes V signature. Assumptions, 
li mitations and numerical t ests of this method are presented 
bv lKochukhov et all l|2010l ). 

We adopted lists of metallic spectral lines suitable 
for earl y-B stars from the Vienna Atomic Line Database 
fVALD. lKupka et al.ir200Gl '). The depth of each line was ad- 
justed to match the observed spectra of each star. In total, 
around 300 lines and 260 lines were used for HD 66665 and 
HD 63425 respectively. All LSD proflles were computed on 
a spectral grid with a velocity bin of 2.6 km s"'^. 

In order to assess the presence of a magnetic signature 
in our LSD proflles, we computed the probability that 
the observed deviation of the signal, with respect toV = G, 
is only due to noise (the false alarm probability, or FAP). 
We consider a signal to be detect ed when this probab ility 
is less than 10~^, as prescribed bv lDonati erakl ([l993). As 
shown in column 5 of Table [T] all observations led to the 
detection of a magnetic signal. The same analysis was per- 
formed on the diagnostic null (A'') profiles, where no signals 
were detected. 

The global longitudinal component of the magnetic field 
was inferred from the first moment of each Stokes V LSD 



profile, as described by I Wade et all l|2000l ). The integration 
range on each side of the Stokes I centre of gravity was set to 
±22kms"^ and ±25kms"^ for HD 66665 and HD 63425 re- 
spectively. While the longitudinal field provides a useful sta- 
tistical measure of the line-of-sight component of the field, 
we stress that we do not use it as the primary diagnostic 
of the presence of a magnetic field. This is because a large 
variety of magnetic configurations can produce a null lon- 
gitudinal field while still generating a detecta ble Stokes V 
profile in the velocity-resolved line profiles (see |Peti"3l201ll ). 

Figure [3] (left panel, right frame) shows the LSD profiles 
of HD 66665. The Stokes V profiles vary on a timescale of 
at most a few days, although a period determination is not 
possible at this point. As the exact rotation phases of our 
ob servations a r e not known, we used the method described 
bv lPetit et al] (|2008| ). which compares the observed Stokes 
V profiles to a large grid of synthetic profiles, described by 
a dipole oblique rotator model. The model is parametrized 
by the dipole field strength Bp, the rotation axis inclination 
i with respect to the line of sight, the positive magnetic 
axis obliquity /? and the rotational phase tp. Assuming that 
only ip may change between different observations of a given 
star, the goodness-of-fit of a given rotation-independent {Bp, 
i, P) magnetic configuration can be computed to determine 
configurations that provides good posterior probabilities for 
all the observed Stokes V profiles, in a Bayesian statistic 
framework. 

In Figure [3] (left panel, right frame) , we show the syn- 
thetic profiles (in solid red) produced by the {Bp, i, /?)- 
configuration that produces the peak of the joint posterior 
probability. The inclined dipole model can reproduce the 
current observations in a satisfactory fashion, although the 
Xrod remains high for two observations (those with the high- 
est s/n). Whether these deviations from the model come 
from artefacts of the LSD procedure, systematics introduced 
by our simple line profile modelling, or a more complex field, 
remains to be investigated once phase-resolved observations 
are obtained and a more detailed analysis can be performed. 

By treating any features that cannot be explained by 
the inclined dipole model as additional Gaussian noise, we 
can obtain a conservative estimate of the model parameters. 
For HD 66665, the posterior probability density marginal- 
ized for Bp (Fig. [3] left panel, left frame) peaks at 670 G, 
and provides a strong constraint on the dipole field lower 
limits (606 G at 99.7%). However, the PDF extends to high 
Bp (1.6 kG at 99.7%), as the exact field strength depends 
strongly on the geometry, as illustrated by the 2D PDFs, 
marginalized for the Bp-i and Bp-j3 planes. 

The LSD profiles obtained for HD 63425 are shown in 
Figure |3] (right panel, right frame). No significant variations 
are observed between the profiles. The red curves show the 
synthetic profiles produced by the {Bp, i, (3) configuration 
that produces the peak of the joint posterior probability. The 
inclined dipole model can once again reproduce the current 
observations in a satisfactory fashion, although some details 
are not reproduced, as shown by the x?cd values. Note that 
our simple modelling of the Stokes I profile was not able 
to reproduce the broad wings of the intensity line profile. 
The posterior probability density marginalized for Bp (right 
panel, left frame) peaks at 460 G, and the 99.7% credible 
region extend from 406 G to 1.8 kG, depending on the exact 
geometry. As no assumptions are made about the phases 
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Figure 3. Magnetic modelling of HD 66665 (left) and HD 63425 (right). For each star, the left panel shows the posterior probability 
density marginalized for the dipole strength Bp (bottom) and for the 2D Bp-i (solid red) and Bp-/3 (dashed blue) planes (top). The 
credible regions containing 68.3%, 95.4%, 99% and 99.7% of the probability are shaded, from dark to light colours respectively. The right 
panel shows the LSD Stokes I profile for one observation (bottom) and all LSD Stokes V profiles (top). The synthetic profiles produced 
by the (Bp, i, /3)-configuration correspond to the peak of the joint posterior probability density, the best phase being extracted for each 
observation. The corresponding X^ed indicated. 



or each observation, the model prefers a dipole field that is 
roughly aligned with the rotation axis (/? ~ 0°), as this con- 
figuration is more likely to occur with observations taken at 
random rotational phases. However, given the low v sini and 
the short time span of our observations, it is possible that all 
our observations were taken at roughly the same rotational 
phase. This would suggest that the stellar rotation period 
is much longer than the period of the observations. Given 
the resemblance of this star with r Sco, and the observed 
variability of the UV line profiles, this hypothesis is likely. 
More observations shall therefore be necessary to constrain 
the exact geometry of the field. 



5 CONCLUSION 

We have presented the characteristics of two stars - 
HD 66665 and HD 63425 - which we believe are analogues 
to the magnetic massive star r Sco for two main reasons: (i) 
The UV spectra of these stars are similar to the once- unique 
spectrum of r Sco. (ii) We have shown that these two stars 
host magnetic fields. 

We have shown that all three of these stars have similar 
fundamental properties. However, the mass-loss rates we es- 
timate for HD 66665 and HD 63425 are lo wer than the value 
adopted for r Sco bv lDonati et~all (|2006l . 6 x 10"* Mq yr'^) 
in their analysis of t hat star. However, that mass-loss rate 
l|Mokiem et al.ll2005h is determined from the optical spec- 
trum of r Sco, and could differ systematically from that de- 
termined from UV line profiles - in particular in the presence 
of a magnetic field. 

The derived mass-loss rates for HD 66665 and HD 63425 
are sufficiently low that they should be consi dered a s stars 
presenting the weak wind problem ((Martins et al.l l2005l : 
iMarcolino et al.ll2009h . since the expected mass-loss rate is 
about 10~^ M0 yr~^ - about 20 time la rger than the infer red 
rates (following the recipe provided bv lVink et al]|2000l ). 

Our modelling of the LSD Stokes V profiles assuming 



an inclined dipole model results in surface field strength^ 
that are comparable to, or perharps slightly larger than, the 
mean surface field strength of r Sco. The current observa- 
tions can be acceptably reproduced by the dipole model, 
although more phase-resolved observations are required in 
order to assess the potential complexity of their magnetic 
fields, and verify if the wind anomalies are linked to the 
field complexity. 
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